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1.  Purpose 

This  engineer  technical  letter  (ETL)  provides  an 
example  of  the  use  of  nonlinear  incremental  structural 
analysis  (NISA)  for  a  massive  concrete  structure 
(MCS). 

2.  Applicability 

This  ETL  applies  to  all  HQUSACE  elements,  major 
subordinate  commands,  districts,  laboratories,  and 
field  operating  activities  (EOA)  having  responsibilities 
for  the  design  of  civil  works  projects. 

3.  References 

a.  EM  1110-2-2006,  Roller  Compacted 
Concrete. 

b.  EM  1110-2-2200,  Gravity  Dam  Design. 

c.  ETL  1 1 10-2-365,  Nonlinear  Incremental 
Structural  Analysis  of  Massive  Concrete  Structures. 

4.  Discussion 

a.  Background.  NISA  is  a  developing  technol¬ 
ogy  which  has  been  used  on  several  recent  civil 
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works  projects  and  valuable  experience  is  gained  with 
each  new  usage  of  the  method.  NISA  is  used  primar¬ 
ily  as  a  supplemental  tool  for  the  design  of  MCS  in 
order  to:  improve  the  cracking  performance  of  a 
particular  type  of  structure  which  has  exhibited  unac¬ 
ceptable  cracking  in  the  past;  to  more  accurately 
predict  the  structural  behavior  of  an  unprecedented 
structure,  or  to  develop  more  cost  effective  structures 
by  revising  the  geometric  configuration,  materials  or 
construction  parameters. 

b.  Application  to  Zintel  Canyon  Dam.  Appen¬ 
dix  A  presents  a  study  conducted  of  the  Zintel  Can¬ 
yon  Dam  in  the  Walla  Walla  District  using  the  NISA 
method  of  analysis.  The  dam  was  constructed  using 
RCC  and  the  NISA  study  was  conducted  by  CENPW 
after  completion  of  construction.  The  study  provides 
insight  into  the  use  of  NISA  as  an  analytical  tool  for 
designers. 

5.  Action 

The  enclosed  study  demonstrates  the  procedures 
necessary  to  conduct  a  NISA  analysis  of  a  RCC 
structure.  NISA  should  be  performed  according  to 
the  requirements  of  ETL  1110-2-365. 


PAUL  D.  BARBER,  P.E. 
Chief,  Engineering  Division 
Directorate  of  Civil  Works 
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APPENDIX  A:  NONLINEAR  INCREMENTAL  STRUCTURAL 
ANALYSIS  OF  ZINTEL  CANYON  DAM 


A-1.  Introduction 

a.  Purpose. 

(1)  The  U.S.  Army  Corps  of  Engineers,  Civil 
Works,  has  been  developing  guidelines  and  proce¬ 
dures  for  determining  thermally  generated  stresses 
resulting  in  construction  of  massive  concrete  struc¬ 
tures.  The  original  product  of  the  research  and  devel¬ 
opment  was  documented  in  ETL  1110-2-324,  Special 
Design  Provisions  for  Massive  Concrete  Structures 
and  has  since  been  superseded  by  ETL  1110-2-365, 
Nonlinear  Incremental  Structural  Analysis  of  Massive 
Concrete  Structures. 

(2)  The  Walla  Walla  District  was  commissioned 
to  perform  a  thermal  stress  analysis  of  Zintel  Canyon 
Dam.  The  goals  in  performing  this  study  are  to: 

(a)  Perform  a  nonlinear,  incremental,  structural 
analysis  (NISA)  to  evaluate  the  effects  of  tempera¬ 
ture,  Roller  Compacted  Concrete  (RCC)  material 
properties,  and  the  subsequent  volume  changes  on  the 
cracking  potential  of  Zintel  Canyon  Dam.  The  pur¬ 
pose  of  such  an  analysis  is  to  evaluate  costs  and 
performance  so  that  appropriate  design  features  and 
requirements  may  be  established.  Obviously,  since 
the  dam  is  complete,  no  design  modifications  will  be 
done.  However,  some  observations  of  the  effective¬ 
ness  of  such  an  evaluation  can  be  made. 

(b)  Demonstrate  the  implementation  of  the  NISA 
process  for  an  RCC  structure  and  compare  analytical 
performance  with  observed  performance  of  the  struc¬ 
ture.  In  the  process,  evaluate  the  NISA  method  to 
provide  recommendations  on  what  measures  may  be 
implemented  to  make  the  process  more  serviceable. 
By  performing  a  NISA  of  a  current  project  at  the 
District  level,  valuable  insight  would  be  gained  on 
whether  this  analytical  method  is  a  suitable  tool  for 
designers  and  implementable  at  the  district  level. 

(c)  Evaluate  the  suitability  of  the  Corps  of  Engi¬ 
neers’  guidance  in  performing  nonlinear,  incremental 
structural  analyses  for  Zintel  Canyon  Dam.  At  the 
time  the  NISA  was  performed,  the  guidance  was 
contained  in  ETL  1110-2-324.  A  new  document, 

ETL  1 1 10-2-265,  has  been  published,  but  most  of  the 
procedures  for  performing  a  NISA  remain  the  same. 


b.  Scope.  This  work  is  limited  to  the  NISA 
evaluation  of  Zintel  Canyon  Dam.  The  analyses 
performed  for  the  project  are  not  intended  to  be  com¬ 
prehensive  evaluations  as  described  in  the  ETL,  but 
abbreviated  evaluations,  more  appropriate  for  a  struc¬ 
ture  of  this  type  and  function.  There  are  two  reasons 
for  this  abbreviation.  The  structure  is  a  simple 
gravity  design,  containing  no  contraction  joints  and 
impounding  no  permanent  reservoir.  Since  transverse 
cracking  of  the  structure  poses  no  threat  to  the  safety 
of  the  structure  or  to  the  public  and  does  not  impact 
the  function  of  the  project,  a  relatively  simple  analy¬ 
sis  is  sufficient.  Secondly,  the  level  of  funding  for 
this  study  was  not  sufficient  to  perform  extensive 
analyses.  Eunding  only  permitted  simple  modeling 
and  limited  evaluations.  However,  the  extent  of 
analysis  was  sufficient  to  evaluate  the  thermal  stress 
performance,  to  provide  recommendations  on  district 
implementation  of  the  NISA,  and  to  evaluate  the 
guidelines  specified  in  the  ETL. 

c.  Report.  This  report  is  organized  to  specifi¬ 
cally  address  the  study  goals.  Section  2  provides 
background  on  the  development  of  finite  element 
analyses  for  mass  concrete  structures  which  led  to  the 
development  of  the  current  procedures  for  performing 
thermal  stress  evaluations  and  Corps  guidance  on  the 
subject.  The  NISA  evaluation  for  Zintel  Canyon 
Dam,  which  includes  the  details  of  performing  the 
evaluation  and  the  project  specific  results,  conclu¬ 
sions,  and  recommendations,  is  contained  in  Section  3 
and  the  referenced  appendices.  Recommendations  are 
provided,  in  Section  4,  on  implementing  a  NISA 
evaluation.  Einally,  recommendations  for  modifica¬ 
tion  of  the  guidance,  to  make  it  more  useful  in  appli¬ 
cation  of  the  NISA  process,  are  contained  in 
Section  5. 


A-2.  Background 

a.  Development  of  the  NISA  process. 

(1)  Mass  concrete  structures  are  different  from 
many  concrete  structures  in  that  material  properties 
have  a  significant  effect  on  the  state  of  stress  in  the 
structure.  These  material  properties  are  not  only 
dependent  on  the  type  and  quantity  of  material,  but 
on  the  age  of  the  concrete,  temperature  of  the 
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concrete,  and  the  state  of  stress  of  the  material. 
Further,  certain  material  properties  exhibit  nonlinear 
performance  and  somewhat  unusual  behavior  at  early 
ages.  Consequently,  definitive  analyses  of  mass 
concrete  structures  require  a  very  complex  analysis 
procedure  requiring  the  definition  of  many  variables. 

(2)  A  further  complicating  factor  is  that  most 
mass  concrete  structures,  such  as  dams,  are  con¬ 
structed  over  a  long  period  of  time.  Consequently, 
thermal  stresses  develop  during  the  construction  phase 
and  may  be  significantly  affected  by  subsequent 
construction  activities. 

(3)  The  evolution  of  an  analysis  package  to 
adequately  model  these  variables  has  been  a  long  and 
tedious  process.  One  of  the  earliest  attempts  to 
develop  an  analysis  system  was  in  1966  during  the 
design  of  Dworshak  Dam.  A  finite  element  analysis 
package  was  developed,  under  contract,  by  research¬ 
ers  at  the  University  of  California,  Berkeley.  This 
system  was  likely  one  of  the  first  such  analysis  tools 
ever  developed.  Over  the  years,  the  finite  element 
analytical  process  has  progressed  to  the  point  where 
many  commercial  vendors  have  provided  quite 
sophisticated  analytical  tools  for  the  evaluation  of 
stress  and  strain  in  concrete  structures.  Unfortu¬ 
nately,  most  of  these  systems  are  general  purpose 
computer  codes  that  do  not  address  the  specific  issues 
of  the  time-dependent  behavior  of  mass  concrete. 

(4)  More  recently,  the  Corps  of  Engineers  has 
initiated  the  development  of  supplemental  codes  and 
techniques  to  enhance  and  refine  the  analytical  pro¬ 
cess.  The  Corps  selected  for  general  use,  ABAQUS, 
a  general  purpose  finite  element  code.  ANATECH 
Research  Corporation,  under  contract  to  the  Water¬ 
ways  Experiment  Station  (WES)  developed  software 
subroutines  to  be  used  with  the  ABAQUS  general 
purpose  finite  element  code.  These  subroutines  were 
designed  to  allow  the  user  to  input  accurate,  time 
dependent  and  cracking  material  properties  of  con¬ 
crete  into  the  ABAQUS  model. 

(5)  ETL  1110-2-324  dated  30  March  1990,  “Spe¬ 
cial  Design  Provisions  for  Massive  Concrete  Struc¬ 
tures,”  was  published,  providing  policy  guidance  to 
designers  for  execution  of  a  NISA  for  the  design  of 
mass  concrete  structures. 

b.  Applicability  of  ETL-1 1 10-2-324.  The  guid¬ 
ance  provided  by  this  ETL  has  been  followed  in  the 
analysis  and  design  of  several  structures.  Most 


notable  was  the  analysis  of  Lock  and  Dam  26  near 
St.  Louis,  Missouri,  and  the  current  design  of 
Olmsted  Locks  and  Dam  on  the  Ohio  River.  Several 
concerns  have  been  raised  regarding  the  implementa¬ 
tion  of  the  ETL.  The  analytical  process,  as  outlined 
in  the  ETL,  is  extremely  comprehensive  and  expen¬ 
sive.  Certainly,  not  all  Massive  Concrete  Structures 
require  a  NISA  to  be  performed.  The  results  of  this 
study  can  be  used  to  determine  the  applicability  of 
the  current  ETL  guidance  for  future  projects  that  may 
or  may  not  require  a  NISA  be  performed. 


A-3.  Zintel  Canyon  Dam  Thermal  Cracking 
Evaluation 

a.  General.  The  purpose  of  performing  a  ther¬ 
mal  cracking  evaluation  for  Zintel  Canyon  Dam  was 
to  determine  the  consequent  cracking  of  the  structure 
resulting  from  thermally  generated  volume  changes. 
Since  the  evaluation  was  performed  some  time  after 
the  construction  of  the  dam,  actual  conditions,  such  as 
ambient  temperatures,  RCC  placing  temperatures,  and 
placing  schedules  were  used  in  the  model.  However, 
additional  laboratory  work  to  characterize  RCC  mate¬ 
rials  was  not  done  because  of  limited  funds.  This 
included  creep  properties,  thermal  properties,  and 
tensile  strain  capacity.  Instead,  these  material  proper¬ 
ties  were  estimated.  In  addition,  the  limited  access  to 
computer  resources  necessary  to  perform  the  analyses 
further  limited  the  depth  of  the  investigation. 

b.  Project  description. 

(1)  Zintel  Canyon  Dam  is  a  straight  axis  con¬ 
crete  gravity  structure.  The  length  of  the  structure  is 
520  ft  across  at  the  crest  and  126  ft  above  the  foun¬ 
dation  at  the  deepest  point.  The  structure  is  con¬ 
structed  of  70,600  cu  yd  of  RCC.  The  outflow 
spillway  is  160  wide  with  a  crest  16  ft  below  the  top 
of  the  dam.  The  spillway  flows  are  contained  by 
cast-in-place  concrete  training  walls  anchored  to  the 
RCC  mass  and  RCC  gravity  training  walls  bordering 
the  stilling  basin.  The  dam,  stilling  basin,  and  stilling 
basin  training  walls  are  founded  on  basalt  rock  (see 
appendix  C  for  a  more  complete  project  description). 

(2)  The  project  provides  flood  protection  to  the 
city  of  Kennewick,  Washington.  It  is  located  on 
Zintel  Canyon,  a  19-square  mile  water  course  which 
threatens  the  city  with  winter  snow  melt  and  summer 
thunderstorm  events.  The  water  course  is  otherwise  a 
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dry  streambed.  The  structure  will  impound  the 
100-year  flood  for  no  more  than  20  days.  A  self¬ 
regulating  outlet  provides  reservoir  drawdown  at  a 
controlled  rate.  The  structure  is  designed  to  require 
no  manned  operations  in  the  event  of  a  flood. 

(3)  After  final  excavation  of  the  foundation,  the 
rock  surfaces  were  cleaned  and  covered  with  a  wet- 
mix  shotcrete  and  foundation  concrete.  RCC  was 
placed  in  12-in.  thick  horizontal  layers  on  and  against 
these  foundation  surfaces.  Interfaces  of  the  RCC  and 
the  foundation  concrete  were  bonded  with  a  bedding 
mortar.  Similarly,  the  RCC  lift  joints  were  fully 
bonded  with  the  same  bedding  mortar.  RCC  place¬ 
ment  began  on  6  July  1992  and  was  completed  on 
15  October  1992.  The  126  RCC  lifts  were  placed  in 
approximately  75  placing  days  during  a  100-day 
period.  In  general,  the  process  was  to  deliver  the 
RCC  to  the  dam  on  a  conveyor.  Front-end  loaders 
received  the  RCC  and  transported  it  to  the  desired 
location.  The  RCC  was  spread  with  a  small  dozer 
and  compacted  with  a  vibratory  roller. 

c.  ABAQUS  model  and  data. 

(1)  After  developing  the  mesh  for  a  two-dimen¬ 
sional  transverse  model  through  the  spillway  section, 
it  became  apparent  that  a  three-dimensional  model 
would  become  extremely  large.  Although  the  CRAY 
could  handle  the  computational  analysis  during  this 
study,  the  time  required  to  execute  such  a  model 
would  have  been  extensive,  due  to  the  system  work¬ 
load.  In  addition,  the  geometry  of  the  dam  does  not 
lend  itself  to  easy  input  generation  for  ABAQUS  and 
as  a  result,  extensive  efforts  to  generate  the  model 
would  be  required.  A  three-dimensional  model  would 
increase  the  cost  of  the  study  significantly  beyond 
initial  estimates.  Because  of  these  factors,  the  study 
was  limited  to  two  two-dimensional  models.  Thermal 
analyses  of  this  nature  have  most  often  been  done 
using  two-dimensional  models.  A  two-dimensional 
transverse  model  usually  gives  good  analytical  results 
that  can  be  used  to  predict  cracks  that  propagate 
inward  from  the  surface  and  cracks  originating  from 
the  foundation.  However,  recent  projects,  most  nota¬ 
bly  Upper  Stillwater  Dam,  experienced  significant 
cracking.  These  cracks  propagate  vertically  from  the 
foundation  and  are  oriented  perpendicular  to  the 
longitudinal  axis.  The  Bureau  of  Reclamation  per¬ 
formed  a  thermal  analysis  to  predict  the  crack 
potential.  They  utilized  a  transverse,  two-dimensional 
model  and  a  longitudinal  model  of  a  horizontal  plane. 


Because  of  the  magnitude  of  the  cracks  at  the 
foundation/RCC  interface,  it  was  necessary  to  deter¬ 
mine  the  suitability  of  a  longitudinal  model  of  a 
vertical  plane  for  use  in  predicting  crack  potential. 

As  a  result,  two  2-dimensional  models  were  analyzed. 
One  model  computed  stresses  resulting  from  a  trans¬ 
verse  model  and  the  other  computed  stresses  from  a 
longitudinal  model  in  a  vertical  plane.  While  these 
two-dimensional  models  provide  less  accurate  solu¬ 
tions  than  a  3 -dimensional  model,  the  accuracy  is 
appropriate  considering  the  size  of  the  project  and 
available  funding  for  the  study.  The  limitation  of  the 
two-dimensional,  longitudinal  model,  is  that  it  is 
assumed  that  symmetry  exists  on  either  side  of  the 
longitudinal  plane,  which  is  not  the  case.  This 
assumption  may  cause  a  shifting  of  the  thermal  gradi¬ 
ent  from  its  actual  location.  In  the  case  of  the  trans¬ 
verse  model,  the  full  section  is  modeled  which  more 
accurately  models  the  thermal  gradient.  Appendix  A, 
Figure  A-1,  shows  the  location  of  the  assumed  two- 
dimensional  plane  through  the  dam.  Thermal  con¬ 
tours  of  both  models,  shown  in  Appendix  A, 
demonstrate  a  good  correlation  of  peak  temperatures 
between  the  two  models. 

(2)  Because  RCC  exhibits  similar  material  char¬ 
acteristics  as  conventionally  placed  mass  concrete,  the 
analytical  process  for  determining  thermal  gradients 
and  stresses  in  RCC  is  practically  the  same.  The 
analytical  procedures  have  been  well  documented  by 
previous  analyses  and  authors  and  should  follow  the 
general  guidance  established  in  ETL  1110-2-324. 
However,  RCC  construction  generally  occurs  over  a 
relatively  short  time  frame  with  numerous  lift  joints 
(usually  1  to  2  ft  in  height).  Conventionally  placed 
mass  concrete  usually  involves  placements  with  lift 
heights  of  4  to  7  ft  with  5-  to  7-day  restrictions 
placed  on  form  removal.  The  exposure  and  depth  of 
each  lift  in  conventionally  placed  concrete  will  gener¬ 
ally  define  the  limits  of  the  number  of  steps  neces¬ 
sary  to  perform  the  incremental  analysis.  In  contrast, 
continuous  placement  of  RCC  on  some  projects  has 
achieved  four  1-ft  lift  heights  in  24  hr.  In  a  7-day 
period,  RCC  can  achieve  changes  in  elevation  of  25 
to  30  ft,  depending  on  the  project  specifics  and  resul¬ 
tant  production  rates.  Therefore,  before  a  NISA  can 
be  performed  for  RCC,  a  comprehensive  study  of 
production  rates  must  be  completed  in  order  to  select 
time  steps,  the  number  of  steps,  and  element  mesh 
size.  The  results  of  the  production  rates,  in  combina¬ 
tion  with  the  element  size,  and  parametric  studies 
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described  in  the  following  paragraph,  will  define  the 
element  mesh. 

(3)  The  mesh  sizes  for  the  two  models  were 
established  by  the  equation  provided  in  the  ABAQUS 
user  manual  and  restated  in  ETL  1110-2-234.  Results 
from  this  equation  were  found  to  be  very  restrictive 
on  the  element  size.  Instead,  a  simple  parametric 
study  was  performed  to  determine  if  a  larger  element 
size  could  be  used  without  creating  numerical  insta¬ 
bilities  in  the  thermal  model.  The  results  of  the  para¬ 
metric  study  indicated  that  the  maximum  length  of 
the  element  in  the  direction  of  heat  flow  could  be 
48  in.  For  both  models,  48  in.  was  the  maximum 
size  of  element  used  in  any  direction  for  a  6-hr  time 
interval.  A  6-hr  time  interval  was  chosen  based  on 
production  rate  of  RCC  and  because  it  satisfied  the 
maximum  time  interval  required  to  compute  early 
heat  gain  in  the  concrete.  Results  of  the  parametric 
study  are  shown  in  Appendix  B.  Mesh  size  was  then 
correlated  with  the  RCC  placement  schedule  of  the 
dam  and  was  designed  to  capture  heat  gains  at  the 
early  ages  of  construction.  Appendix  A,  Figure  A-1 
represents  a  time  history  for  construction  and  the 
analysis  for  the  transverse  and  longitudinal  model 
respectively,  as  well  as  the  initial  conditions,  bound¬ 
ary  conditions,  and  input  values  used.  Boundary 
conditions  include  the  insulating  effects  of  upstream 
precast  facing  panels,  free  surface  convection,  soil 
(rock)  conditions,  downstream  stilling  basin  slab, 
average  daily  temperatures  for  preconstruction,  during 
construction,  and  postconstruction.  Initial  conditions 
include  RCC  placement  temperatures,  and  initial 
foundation  temperatures.  The  user  subroutine 
DFLUX  was  used,  in  conjunction  with  ABAQUS,  to 
generate  time-dependent  heat  fluxes  for  the  thermal 
analysis.  Parameters  used  in  DFLUX  included  adia¬ 
batic  heat  gain  (time  and  temperatures)  for  the  RCC 
mix  as  well  as  initial  placement  times  for  each  lift. 
Adiabatic  heat  gain  curve  is  plotted  in  Figure  A-1. 
Results  of  the  thermal  analysis  for  both  the  transverse 
and  longitudinal  model  are  represented  in  contour 
plots  in  Appendix  A,  Figures  A-2  to  A- 11  and  time 
history  plots  of  maximum  nodal  temperatures  in 
Appendix  A,  Figures  A- 12  and  A- 13.  The  maximum 
temperature  reported  in  the  transverse  model  is 
represented  by  node  2330,  Figure  A-12c,  and  in  the 
longitudinal  model  by  node  3213,  Figure  A- 13c. 
Maximum  temperatures  and  temperature  differential 
correlate  well  with  predicted  temperatures  calculated 
by  using  approximate  computational  methods  for 
Zintel  Canyon  Dam.  As  stated  previously,  effects  of 
maximum  heat  gain  between  the  two  models  were 


nearly  the  same.  However,  there  is  a  notable  differ¬ 
ence  in  the  rate  of  coding.  The  two-dimensional 
elements  do  not  have  the  capability  to  conduct  heat  in 
the  out-of-plane  direction.  With  convection  being 
modeled  along  the  top  surface  only,  of  the  longitudi¬ 
nal  model,  the  elements  sustain  a  higher  thermal 
gradient  over  a  much  longer  period  of  time. 

(4)  Discussion  of  stresses  for  both  models  should 
be  limited  to  principle  tensile  stresses.  However, 
software  developed  for  plotting  stress  histories  is  only 
capable  of  plotting  stress  in  orthogonal  directions  and 
for  shear  stresses.  Because  Zintel  Canyon  Dam  nor¬ 
mally  has  no  pool  and  experiences  a  very  short 
duration  reservoir  impoundment,  cracking  posed  no 
concerns  related  to  seepage.  Of  concern  are  the 
orientation  of  cracks  that  will  compromise  the  stabil¬ 
ity  of  the  structure.  A  feature  of  the  ABAQUS-based 
NISA  that  sets  it  apart  from  others  is  that  it  allows 
material  properties  and  relationships  to  be  user- 
defined.  UMAT  is  the  subroutine  that  provides  a 
time-dependent  cracking  material  model.  The  subrou¬ 
tine  allows  input  of  specific  material  properties  and 
calibration  of  the  model-predicted  properties  against 
actual  observed  material  performance.  While  some 
material  testing  had  been  done  for  Zintel  Canyon 
Dam,  extensive  evaluation  of  time-dependent  proper¬ 
ties  and  creep  performance  had  not  been  determined. 
Consequently,  calibration  of  the  UMAT  material 
model  could  not  be  performed.  The  material  model 
generated  for  Olmsted  Locks  and  Dam  was  used 
except  that  the  material  constants  were  replaced  with 
actual  or  estimated  values  for  Zintel  Canyon  Dam 
RCC  materials.  This  means  that  the  UMAT  predicted 
performance  for  Zintel  Canyon  Dam  was  based  on 
Olmsted  material  relationships.  No  data  were  avail¬ 
able  to  shift  the  relationship  curves.  Without  actual 
data  to  calibrate  the  material  model,  changes  would 
be  arbitrary  and  not  necessarily  an  improvement  over 
using  the  Olmsted  data.  To  perform  a  reliable  NISA, 
these  material  properties  need  better  definition.  More 
complete  definitions  of  the  development  of  the  modu¬ 
lus  of  elasticity  and  creep  are  critical  to  accurate 
results.  For  smaller  scale  projects,  standard  relation¬ 
ships  for  a  range  of  materials  should  be  developed  so 
that  the  analyst  can  select  the  performance  relation¬ 
ship  that  most  likely  models  the  materials  being  eval¬ 
uated.  Only  the  larger  projects  will  have  the  funding 
to  perform  a  complete  battery  of  laboratory  evalua¬ 
tions.  For  the  longitudinal  model,  the  maximum 
principle  stress  occurred  at  the  foundation/RCC  inter¬ 
face  at  element  1796  and  are  presented  by  stress 
contour  plots  in  Appendix  A,  Figures  A- 14  to  A- 19 
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and  A-20  to  A-25.  ABAQUS  calculated  a  maximum 
principle  stress  of  715  psi  which  occurs  during  the 
cooling  period  when  maximum  temperature  differen¬ 
tials  occur.  Earlier  testing  of  the  RCC  mix  indicated 
a  28-day  tensile  strength  capacity  of  200  psi.  The 
715-psi  principle  stress  calculated  is  far  in  excess  of 
the  tested  direct  tensile  capacity.  This  may  be  attrib¬ 
uted  to  the  fact  that  the  analyses  were  performed 
using  the  original  version  of  UMAT.  That  version 
contained  inconsistencies  which  resulted  in  the  crack¬ 
ing  threshold  to  be  computed  in  an  unconservative 
manner.  This  inconsistency  has  been  corrected  in  the 
current  version  of  UMAT  which  was  not  available 
when  this  study  began.  However,  newer  versions  of 
software  available  at  WES  have  the  capability  to 
predict  and  plot  direction  and  magnitude  of  cracks. 
Observation  of  the  stress  contours  for  the  longitudinal 
model  indicate  higher  stresses  occur  at  the  foundation 
interface  and  in  the  upper  reaches  of  the  abutment 
and  at  the  spillway.  This  can  be  attributed  to  the 
temperature  differential  occurring  between  the 
exterior  and  interior  elements  and  the  rigidity  of  the 
foundation.  At  the  time  of  year  for  postconstruction 
cooling,  the  average  ambient  temperature  is  decreas¬ 
ing  causing  a  larger  temperature  differential  from 
surface  elements  to  interior  elements.  Time  history 
stress  plots  of  various  points  of  high  stresses 
observed  in  the  contour  plots  are  presented  in  Appen¬ 
dix  A,  Eigures  A-26  to  A-32  and  A-33  to  A-40.  The 
stress  time  histories  presented  in  Eigures  A  10.4  and 
A10.6  indicate  that  some  cracking  may  have 
occurred.  Jumps  in  stress,  as  indicated  in  these  plots, 
typically  do  not  occur  unless  a  crack  has  formed. 
Discussion  of  these  results  are  presented  below. 
Likewise,  the  transverse  model  predicted  stresses  that 
are  higher  than  the  limited  cracking  stress.  The  high¬ 
est  stress  occurred  at  element  217  at  the  foundation/ 
RCC  interface.  In  addition,  a  region  of  high  stress 
occurs  along  the  downstream  exposed  face  of  the 
dam. 

d.  Discussion. 

(1)  A  simple  thermal  analysis  was  done  for  the 
project  during  the  design  phase  of  the  project.  This 
analysis  indicated  that  during  the  normal  summer 
weather  conditions  at  the  site,  the  structure  may  crack 
at  three  locations.  Two-crack  locations  were 
estimated  to  be  located  where  the  foundation  changes 
from  a  horizontal  surface  at  elevation  635  to  the 
sloped  abutments.  The  third  crack  was  speculated  to 
be  in  the  center  of  the  spillway.  None  of  these 
cracks  pose  a  threat  to  the  stability  of  the  structure 


since  the  orientation  will  be  in  the  traditional 
upstream-downstream  direction.  In  addition,  since  the 
dam  is  almost  continuously  dry,  and  channel  flows 
carry  a  phenomenal  silt  load,  this  cracking  of  the 
structure  is  not  of  great  concern.  No  additional 
expense  was  warranted  in  lessening  the  cracking 
potential  by  reducing  the  placing  temperatures  or  by 
installing  transverse  joints. 

(2)  Postconstruction  inspections  revealed  one 
crack  in  the  structure  located  high  on  the  left  abut¬ 
ment  as  a  result  of  a  slope  change  in  the  foundation. 
No  other  cracks  are  apparent.  Less  cracking  has  been 
observed  because  the  restraint  provided  by  the  foun¬ 
dation  is  probably  lower  than  full  restraint  assumed 
by  the  simple  analysis. 

(3)  Examination  of  the  temperature  history  plots 
indicates  that  the  longitudinal  model  cools  at  a  much 
slower  rate.  The  benefits  of  performing  a  longitudi¬ 
nal  model  are  significant  when  time  and  costs  are  of 
concern  for  smaller  projects,  or  projects  of  this  type 
where  certain  cracking  will  not  adversely  effect  the 
performance  of  the  structure.  However,  in  this  case, 
the  thermal  gradients  of  the  longitudinal  model  should 
be  calibrated  with  the  more  accurate  transverse  model 
to  produce  nearly  the  same  rate  of  cooling.  This  may 
become  a  significant  factor  in  the  analysis  because: 

•  While  cooling  is  at  a  slower  rate  than  might 
be  expected,  this  will  cause  volume  changes 
and  hence  the  maximum  stresses  to  occur  at 
a  later  age  in  the  model.  Since  the  higher 
stresses  occur  later  in  time,  the  aging  modu¬ 
lus  will  be  higher.  Since  the  criteria  for 
cracking  is  partially  based  on  stress,  the 
potential  for  cracking  will  be  unconservative 
for  this  case. 

•  The  thermal  stresses  are  being  applied  at  a 
later  age;  therefore,  the  effect  of  creep  will 
play  a  less  significant  role  in  stress  relax¬ 
ation.  This  may  be  conservative;  however, 
as  long  as  a  large  amount  of  effort  has  been 
expended  to  accomplish  a  material  investiga¬ 
tion,  it  would  be  prudent  to  spend  the  same 
amount  of  effort  to  calibrate  the  longitudinal 
model.  This  will  ensure  that  the  analysis  will 
incorporate  the  more  accurate  creep  data  in 
the  earlier  time  steps  when  cooling  would  be 
expected  to  cause  higher  stresses.  Hence,  the 
effects  of  stress  relaxation  due  to  creep  will 
be  incorporated  into  the  analysis. 
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(4)  Calibrating  the  longitudinal  model  may  be 
done  by  adjusting  the  thermal  conductivity  of  differ¬ 
ent  element  sets,  allowing  for  higher  conductivity  in 
the  earlier  time  steps  and  reducing  the  conductivity  in 
the  later  time  steps  where  the  convective  surface 
plays  a  more  significant  role  in  cooling. 

(5)  Examination  of  the  stress  history  plots  indi¬ 
cates  stresses  are  still  increasing  as  a  result  of 
decreasing  ambient  air  temperatures.  For  both  mod¬ 
els  the  cooling  period  should,  at  a  minimum,  be 
applied  until  thermal  stresses  begin  to  decrease  and 
attain  a  steady  state.  For  this  analysis  the  cooling 
period  was  9  months  for  the  transverse  model  and 

3  months  for  the  longitudinal  model.  Admittedly,  this 
analysis  fell  short  of  predicting  maximum  stresses. 
However,  review  of  the  principle  stress  contours  from 
the  longitudinal  model  reveals  high  stress  areas  where 
cracking  is  likely  to  occur.  One  area  is  at  the  center 
of  the  spillway,  the  second,  at  the  intersection  of  the 
spillway  and  non-overflow  section,  and  the  third,  in 
the  upper  reaches  of  the  abutments.  With  the  version 
of  UMAX  used  for  this  study,  it  is  difficult  to  predict 
where  cracking  will  initiate  first.  However,  where 
locations  of  high  stresses  occur,  the  model  may  be 
modified  to  depict  the  location  of  transverse  joints  if 
cracking  is  undesirable  in  those  regions.  Other  condi¬ 
tions  that  affect  the  stress  in  the  dam  are  the  assumed 
foundation  restraints.  Both  models  include  fully 
restrained  boundary  conditions  at  the  RCC/foundation 
interface.  By  visual  observations  of  the  rock  and 
postconstruction  coring  of  the  foundation,  the 
assumed  restraint  conditions  used  in  the  model  could 
be  modified  to  provide  a  more  flexible  restraint  con¬ 
dition  or  an  adjusted  foundation  modulus.  Before 
proceeding  with  any  further  analysis,  the  most  recent 
version  of  UMAX  should  be  used  to  include  the 
redistribution  of  stress  that  occurs  after  cracking.  For 
this  study,  no  further  calibration  of  the  longitudinal 
thermal  model  was  completed.  Xhis  is  mainly  due  to 
limited  scope  of  the  study  and  the  fact  that  cracking 
is  not  of  great  concern  for  Zintel  Canyon  Dam. 

e.  Conclusions. 

(1)  Performance  of  the  analysis  leads  to  several 
conclusions  and  recommendations  for  subsequent 
steps  to  proceed  with  further  evaluation  of  cracking. 
Xhe  first  step  would  be  to  make  the  required  adjust¬ 
ments  in  the  foundation  modulus  and/or  the  restraint 
conditions.  Xhe  model  for  Zintel  Canyon  Dam 
included  the  foundation  in  the  stress  analysis. 
Adjusting  the  foundation  modulus  to  tested  values 


would  be  more  appropriate  in  this  case.  If  the  foun¬ 
dation  conditions  are  modeled  by  the  use  of  springs, 
the  degree  of  restraint  may  be  adjusted  by  softening 
or  stiffening  the  spring  constants.  In  the  case  of 
Zintel  Canyon  Dam,  the  foundation  rock  was  highly 
fractured.  Xherefore,  the  degree  of  restraint  provided 
by  the  concrete  to  rock  interface  may  be  less  than 
fully  restrained.  However,  unless  sufficient  data 
support  softening  the  spring  constants,  the  results  of 
the  stress  analysis  would  be  unconservative  in  this 
case.  Secondly,  the  model  can  be  modified  by  insert¬ 
ing  a  transverse  joint  at  the  mid-point  of  the  spillway 
or  at  the  comers  of  the  spillway/non-overflow  inter¬ 
section.  Additional  measures  could  include  reducing 
the  RCC  placement  temperature  and  mandating  place¬ 
ment  schedules  to  avoid  hot  seasons. 

(2)  For  Zintel  Canyon  Dam,  considering  the 
frequency  of  reservoir  impoundment  and  the  function 
of  the  structure,  the  observed  cracking  is  well  within 
acceptable  levels  for  the  project.  Xhere  would  be  no 
benefit  for  implementing  these  measures. 

(3)  Xhe  results  of  this  NISA  analysis  were,  in 
general,  consistent  with  the  results  of  the  approximate 
thermal  analysis  performed  during  design  of  the  pro¬ 
ject.  Both  indicate  that  cracking  may  occur  in  three 
areas.  Based  on  that  comparison  and  the  observed 
performance  of  Zintel  Canyon  Dam  and  other  RCC 
dams,  the  NISA  for  Zintel  Canyon  Dam  provided  no 
additional  information  to  attain  the  desired  objectives 
stated  in  the  EXF.  Xhis  statement  is  based  on  the 
fact  that: 

•  Some  cracking  would  be  acceptable  as  long 
as  structural  stability  is  not  compromised. 

•  Joints  for  control  of  cracking  is  not  necessary 
for  serviceability  conditions  since  there  is  no 
permanent  reservoir. 

•  Because  of  the  combinations  of  the  above, 
real  cost  savings  have  been  achieved  by 
maximizing  production  rates  for  placement  of 
RCC. 

•  Unusual  loadings,  extreme  loadings,  or 
severe  operational  conditions  do  not  exist. 

(4)  Xhis  is  not  to  say  that  a  NISA  should  not  be 
performed  for  RCC  structures.  Each  structure  is 
subject  to  unique  conditions  and  loadings  resulting  in 
unique  structural  features.  Xhese  factors  are 
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evaluated  and  developed  by  a  team  of  responsible 
engineers  who  must  determine  the  level  of  analyses 
and  consultation  necessary  to  obtain  any  of  the 
desired  objectives  stated  in  the  ETL.  The  ETL  can¬ 
not  provide  guidance  for  all  possibilities. 

A-4.  NISA  Application  Recommendations 

a.  Applications. 

(1)  This  work  would  not  have  progressed  with¬ 
out  the  aid  of  several  WES  personnel.  These  individ¬ 
uals,  currently  running  ABAQUS  problems  on  the 
CRAY,  provided  a  significant  amount  of  consultation. 
This  type  of  aid  will  likely  be  necessary  should  dis¬ 
trict  users  begin  to  utilize  this  software  and 
equipment. 

(2)  This  evaluation  was  performed  over  a 
lengthy  period  of  time.  The  time  period  was  much 
longer  due  to  a  start-stop  approach  employed  for  the 
study.  Several  problems  were  identified  during  this 
period  that  seem  to  be  shortcomings  of  the  current 
NISA  process.  Admittedly,  some  of  the  problems 
encountered  were  due  to  the  protracted  approach 
taken.  A  start-stop  operation  is  rarely  an  efficient 
operation. 

(3)  Accessing  the  CRAY  computer  provides  a 
variety  of  problems.  Remote  access  for  field  use  of 
ABAQUS  at  this  time  is  not  possible,  because 
ABAQUS  is  currently  site  licensed  for  use  at  WES. 
Districts  preparing  to  embark  on  such  a  study  will 
find  it  necessary  to  utilize  WES  personnel  at  WES  or 
perhaps  negotiate  a  contract  with  the  ABAQUS 
owner  (Hibbitt,  Karlsson,  and  Sorenson  Inc.)  to  allow 
the  District  access  to  ABAQUS  via  the  WES  Cray. 
Our  conclusion  is,  that  while  the  CRAY  may  provide 
significant  computing  capability,  the  logistical  prob¬ 
lems  to  off-site  users  may  provide  more  problems 
than  solutions  for  off-site  users. 

(4)  There  has  been  a  tremendous  amount  of 
effort  expended  in  developing  NISA  via  ABAQUS 
and  the  user  subroutines.  The  logical  next  step  is  to 
develop  desktop  software  packages  that  would  be 
beneficial  to  Districts  with  smaller  projects  requiring 
finite  element  analysis  or  larger  projects  where 
preliminary  two-dimensional  modeling  is  necessary 


prior  to  embarking  on  a  full  scale  three-dimensional 
NISA.  The  Cray-based  software  usage  would  then  be 
reserved  for  those  rare  monumental  projects.  It  is 
recommended  that  a  micro-based  version  of  this 
software  be  utilized  for  most  applications. 

(5)  A  critical  deficiency  in  the  ABAQUS  analyti¬ 
cal  package  is  the  lack  of  graphic  preprocessing  and 
postprocessing  routines.  Adequate  preprocessing 
would  eliminate  much  of  the  input  generation  errors. 
Currently,  there  is  no  interactive  preprocessing  capa¬ 
bility.  Eurther  development  of  preprocessing  and 
postprocessing  routines  should  match  formats 
provided  by  some  of  the  more  common  processing 
routines  currently  being  used  by  designers. 

(6)  Eor  these  analyses,  input  to  the  user  subrou¬ 
tines  required  changing  the  EORTRAN  code  to  pro¬ 
vide  heat  flux  information,  creep  and  shrinkage 
characteristics  corresponding  to  appropriate  element 
sets.  Although  many  engineers  can  decipher 
EORTRAN  code,  many  involved  in  the  work  may  not 
have  produced  any  programs  for  years.  This  can  be 
time  consuming  and  is  a  waste  of  effort  in  the  design. 
A  more  appropriate  means  for  entering  data  would  be 
via  batch  files  to  the  user  subroutines  as  is  currently 
permitted.  UMAT  is  an  extremely  comprehensive 
subroutine  that  is  the  crux  of  the  time-dependent 
stress  analysis.  Subsequent  to  this  study,  the  capabil¬ 
ity  for  entering  data  to  the  user  subroutine  through  a 
batch  file  was  implemented. 

(7)  Eor  Districts  to  utilize  ABAQUS  in  order  to 
execute  a  NISA  as  part  of  the  design  process,  and  to 
achieve  the  stated  objectives  in  the  ETL,  ABAQUS 
must  be  readily  available  to  District  designers.  With¬ 
out  these  labor  saving  additions,  ABAQUS  use  by 
designers  outside  of  WES  may  never  develop. 

(8)  Wide  usage  by  Districts  will  no  doubt  require 
sophisticated  support  services  in  the  form  of  training 
and  user  support.  This  service  is  best  provided  by 
WES.  Program  orientation  is  recommended  in  the 
form  of  a  periodic  PROSPECT  course  on  the  use  of 
the  NISA/ABAQUS  system.  Secondly,  a  staff  mem- 
ber(s)  must  be  available  to  service  the  ABAQUS 
system  and  provide  user  support.  User  support  may 
range  from  troubleshooting  user  problems  to  working 
as  part  of  the  design  team. 
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A-5.  ETL  Recommendations 

a.  Future  guidance. 

(1)  Future  guidance  concerning  mass  concrete 
should  address  whether  performing  a  NISA  is  manda¬ 
tory.  Current  guidance  is  unclear  i.e.,  a  reader  can 
either  assume  that  a  NISA  is  mandatory  (ETL  1110- 
2-324  paragraph  la)  or  that  the  need  for  such  an 
investigation  is  subject  to  consideration  (ETL  1110-2- 
324  paragraph  lb).  It  should  be  noted  that  a  new 
ETL  has  been  published,  that  incorporates  updated 
information  based  on  NISA’s  that  have  been  per¬ 
formed  by  WES  that  may  address  these  issues.  The 
considerations  for  when  to  implement  a  NISA  versus 
other  less  comprehensive  analyses  need  to  be  devel¬ 
oped  and  included  in  any  comprehensive  policy 
document. 

(2)  Certain  basic  questions  must  be  addressed 
prior  to  embarking  on  a  NISA.  These  include: 

•  Why  do  a  thermal  study? 

•  Is  a  thermal  study  appropriate  for  this  struc¬ 
ture?  If  so,  what  level  of  analysis  is 
necessary? 


•  What  are  the  basic  principles? 

•  How  do  I  do  a  thermal  study? 

(3)  It  is  recommended  that  future  documentation 
address  these  issues.  The  document  should  address 
the  basic  issues  of  the  goals  and  desired  objectives 
when  performing  thermal  studies  and  supplemented 
with  more  specific  information,  for  all  levels  of 
analyses. 

(4)  Guidance  in  the  ETL  references  acceptable 
bandwidths  to  be  applied  to  the  mechanical  properties 
of  the  concrete  for  estimated  data.  There  are  several 
ways  to  generate  reasonable  estimates  of  this  data  at 
the  time  of  the  analysis.  Eurther  guidance,  that  would 
be  beneficial,  should  reference  other  sources  that 
contain  the  methodologies  to  estimate  the  data. 

(5)  We  concur  with  the  recommendation  of  the 
ETL  that  the  full  intended  benefit  of  a  NISA  requires 
the  combined  efforts  of  the  structural  designer,  mate¬ 
rials  engineer,  cost  engineer,  and  the  geotechnical 
engineer. 
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Figure  A-1.  Finite  eiement  mesh  and  modei  data 
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Figure  A-4.  Temperature  contours  transverse  model,  ZIntel  Canyon  Dam  12  placement 
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ZINTEL  CANYON  DAM  12  PLACEMENT 
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Figure  A-10.  Temperature  contoura  longitudinal  model,  Zintel  Canyon  Darmtop  ptacement 
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Figure  A-12.  Transverse  temperature  hisiortos 
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A-13.  Noda)  tflitipvralure  histories 
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Figure  A-16.  Principle  Stress  Contours  Transverse  Model,  Zintel  Canyon  Dam,  cool-down  period 
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Figura  A.17.  {global  “X"  Stress  Contours  Trsntvsrss  Model,  Ziniel  Canyon  Da<n,  Dam  20  placement 
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Figure  A-IB.  Global  "Y"  Strees  Contour*  Transverse  Model,  Ztnlel  Canyon  Dam,  Dam  20  pjacement 
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Figure  A-19.  Principle  Stress  Contours  Transverse  Modal,  Zlntal  Canyon  Dam,  Dam  20  placement 
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A-20.  Global  "X”  Stress  Contoura  Longitudinal  Model,  ZInte)  Canyon  Dam,  dam-top  placement 
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Figurs  A-21.  Globa J  "Y*  Strassea  Longitudinal  Mkxlel,  ZInlal  Canyon  Dam,  dam-top  placament 
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a.  Stress  history  element  4B9,  global  X  stresses,  transverse  analysis 


b.  Stress  history  element  489,  global  Y  stresses,  transverse  model 


Time  (Days) 


c.  Stress  history  element  489,  shear  stresses,  transverse  analysis 


Time  (Days) 


Figure  A-31.  Stress  Histories  for  Element  489 
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a.  Stress  history  element  92.  global  X  stresses,  longitudinal  model 


b.  Stress  history  element  92,  global  Y  stresses,  longitudinal  modd 


Time  (Dayel 


c.  Stress  history  element  92,  shear  stresses,  longitudinal  model 


Figure  A-33.  Stress  Histories  for  Element  92 
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a.  Stress  history  dement  518,  global  X  stresses,  longitudinai  mtxlel 


b.  Stress  history  dement  518,  gtobal  Y  stresses,  longitudinal  model 


c.  Stress  history  element  518,  shear  stresses,  longitudinal  model 


Figure  A-35.  Stress  Histories  for  Element  518 
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APPENDIX  B:  ELEMENT  SIZE 


B-1 .  Parametric  Study 

a.  General.  The  integration  procedures  used  in 
the  program  for  transient  heat  transfer  analysis 
requires  a  relationship  between  the  minimum  time 
step  and  the  element  size.  The  equation  to  establish 
this  relationship  is  given  as: 

Af  >  (pc/6A:)A/^  or  <  I  J(6A:At/pc)  (B-1) 

where: 

At  =  time  step 
p  =  density 
c  =  specific  heat 
k  =  thermal  conductivity 
Al  =  element  dimension 

b.  Adiabatic  heat  gain.  Adiabatic  heat  gain  in 
concrete  begins  within  the  first  12  hr  after  placement 
and  can  continue  rapidly  until  a  maximum  is  attained. 
Therefore,  when  performing  incremental  time  depen¬ 
dent  stress  analysis  for  concrete,  it  is  important  to 
keep  the  time  steps  sufficiently  small  during  the  early 
stages  of  the  analysis.  Input  of  the  appropriate  prop¬ 
erties  for  Zintel  Canyon  Dam  into  Equation  B-1 
yields  a  maximum  length  of  element,  using  a  6-hr 
time  interval,  of  27  in.  Analysis  for  a  12-hr  time 
interval  yields  a  38-in.  element.  A  12-hr  time  inter¬ 
val  is  not  a  good  choice  for  calculating  early  heat 
gain  in  the  concrete,  while  placing  a  27 -in.  restriction 


for  a  6-hr  time  interval  doubled  the  size  of  the  model. 
A  48-in.  step  height  nearly  matched  production  rates 
for  daily  concrete  placements,  however,  did  not  fit  the 
criteria  established  in  ABAQUS.  Therefore,  this 
study  focused  on  a  48-in.  element  size  to  determine 
its  reliability  in  reporting  temperature  data. 

c.  One-dimensional  heat  flow.  This  study  is  a 
simple  one-dimensional  heat  flow  problem,  using 
material  properties  for  Zintel  Canyon  Dam.  Two 
models  were  generated,  one  with  a  24-in.  element 
size  in  either  direction  and  the  other  with  a  48-in. 
element  size  in  either  direction.  Depicted  in  Fig¬ 
ure  B-1  are  the  two  finite  element  meshes,  and  boun¬ 
dary  conditions  used  for  the  study.  One  exterior  row 
of  boundary  nodes  were  held  at  a  constant  50  deg 
while  the  ambient  surface  conditions  along  the  oppo¬ 
site  face  was  a  fixed  90  deg.  The  thermal  models 
calculated  nodal  temperatures  in  0.25  day  increments 
for  a  period  of  10  days.  Corresponding  nodal  tem¬ 
peratures  were  compared  from  both  models  to  deter¬ 
mine  accuracy,  and  if  stable  heat  gain  was  being 
generated.  Figure  B-2  contains  plots  of  nodal  tem¬ 
peratures  for  both  the  24-  and  48-in.  meshes  for 
various  times.  The  only  inconsistency  was  at  time  t 
=  0.5  days,  for  the  48-in.  mesh,  where  a  slight  incon¬ 
sistency  in  the  heat  gain  exists.  This  can  be  seen  in 
Figure  B-2.1.  Figures  B-2a  through  B-2d,  indicate 
nearly  identical  heat  gain,  when  comparing  nodal 
temperatures  at  the  same  time  steps  of  the  two 
models.  Because  this  amount  of  inconsistency  was  at 
small,  and  only  occurred  at  one  time  step,  it  was 
considered  negligible  and  would  not  effect  the  out¬ 
come  of  the  study.  Therefore,  we  decided  to  use  a 
time  step  of  6  hr  (0.25  days)  and  a  maximum  element 
size  in  any  direction  of  48  in. 
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Figure  B-1.  Finite  element  mesh  and  input  data 


Temperature  plots  for  t  -  0.6  days 
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Figure  B-2.  Temperature  comparisons 
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APPENDIX  C:  ZINTEL  CANYON  PROJECT 


C-1 .  Project  Description 

a.  General.  The  Zintel  Canyon  Project 
(Figures  C-1  through  C-4)  was  authorized  for  con¬ 
struction  by  resolution  of  the  House  and  Senate  Com¬ 
mittees  on  Public  Works,  December  1970,  under 
authority  of  section  201,  Flood  Control  Act  of  1965 
(Public  Law  298,  89  Congress).  The  project  was 
constructed  substantially  as  authorized.  Detention 
storage  was  reduced  from  2,560  to  1,260  acre-ft  since 
this  was  considered  the  optimum  economic  size  of  the 
dam.  This  alternative  will  not  prevent  damages  in 
some  areas  of  Kennewick  or  avoid  the  use  of  streets 
as  a  channel  during  flooding  in  excess  of  50  years 
(100-year  thunderstorm).  Zintel  Canyon  Project 
includes  a  dam  and  a  floodway  channel  with  required 
structures  that  carry  the  combined  flows  from  the 
dam  and  areas  below  the  dam  through  a  developed 
section  of  Kennewick  to  a  discharge  point  at  the 
Columbia  River.  Zintel  Canyon  Dam  is  a  90-ft  reten¬ 
tion  straight  axis  concrete  gravity  structure  totaling 
approximately  70,000  cu  yd  of  roller  compacted 
concrete  (RCC).  The  purpose  of  the  dam  is  to  pro¬ 
vide  flood  protection  to  the  city  of  Kennewick,  Wash¬ 
ington,  by  impounding  flood  flows  behind  the  dam  up 
to  the  100-year  return  frequency,  and  discharging  that 
volume  over  a  20-day  period.  The  Floodway  Chan¬ 
nel  improvements  consists  of  a  buried  conduit 
designed  to  pass  up  to  a  50-year  composite  flood 
level.  The  78-in.  buried  conduit  is  designed  to  carry 
400  cfs  from  its  intake  at  West  7th  Avenue  and  Van¬ 
couver  to  the  outlet  in  the  Tri-City  Country  Club  Golf 
Course.  From  there  the  natural  channel  is  designed 
to  pass  620  cfs  through  the  Burlington  Northern 
railroad  fill  (Figure  A-1).  Downstream  of  the  railroad 
fill  the  channel  is  designed  to  provide  standard  proj¬ 
ect  protection.  The  project  is  co-funded  by  the 
U.S.  Army  Corps  of  Engineers  (75  percent)  and  the 
city  of  Kennewick  (25  percent).  The  project  is 
located  in  a  semi-arid  region  of  eastern  Washington 
and  borders  on  the  south  end  of  Kennewick,  Wash¬ 
ington.  The  basin,  a  well  defined  water  course  called 
Zintel  Canyon,  is  normally  dry  and  drains  approxi¬ 
mately  28  square  miles  of  the  north  side  of  the  Horse 
Heaven  Hills  of  which  approximately  19  square  miles 
in  area  is  upstream  of  the  dam.  The  drainage 
upstream  of  the  dam  collects  winterstorm  and  thun¬ 
derstorm  runoff,  thereby  providing  a  100-year  flood 
storage  volume  of  1,260  acre-ft. 


b.  Geology  and  foundation.  Zintel  Canyon  is 
located  on  the  southwest  flank  of  the  Pasco  Basin,  a 
structural  feature  formed  by  downward  folding  and 
faulting  of  the  Columbia  River  Basalt  formation. 
Erosion  and  deposition  has  modified  the  structural 
features  by  partially  filling  the  basin  with  sediments 
and  covering  the  rock  slope  with  a  mantel  of  fine¬ 
grained  materials.  Bedrock  is  close  to  the  surface 
within  the  drainage  area  of  Zintel  Canyon  and  where 
the  dam  was  located.  The  foundation  rock  was  com¬ 
posed  of  hard,  dense  basalt  with  closely  spaced  frac¬ 
tures.  The  moderately  un weathered  pieces  were 
bounded  by  weathered  fracture  surfaces.  Eracture 
fillings,  particularly  near  the  surface,  were  filled  with 
silt  and  clay.  Because  the  rock  would  easily  dislodge 
when  the  joint  filling  dried,  as  well  as  from  subse¬ 
quent  construction  activities,  the  exposed  foundation 
rock  was  covered  with  a  minimum  8-in.  layer  of 
pumped  concrete  prior  to  RCC  placement. 

c.  Dam,  spillway  and  outlet.  The  dam  is  a 
straight  axis  concrete  gravity  structure  with  a  crest 
length  of  520  ft  and  a  160-ft,  centrally  located, 
ungated  overflow  spillway.  The  height  of  above 
foundation  is  126  ft  and  86  ft  above  the  existing 
channel  with  a  20-ft  crest  width  in  the  nonoverflow 
section.  The  slope  of  the  downstream  face  was 

.85  horizontal  to  1  vertical  to  facilitate  free  forming 
of  the  downstream  face.  The  upper  24  ft  of  the 
downstream  face  of  the  dam  (adjacent  to  the  spill¬ 
way)  was  constructed  using  vertical  concrete  facing 
panels  as  was  the  upstream  face.  An  80-ft  long 
hydraulic  jump-type  stilling  basin  was  located  at  the 
toe  of  the  structure.  This  stilling  basin  consists  of  a 
12-ft-thick  RCC  base  slab  integrally  constructed  with 
an  RCC  endsill  and  RCC  gravity  training  walls.  The 
spillway  was  designed  to  discharge  a  flow  of 
38,950  cfs.  The  full  width  of  the  spillway  crest  was 
surfaced  with  a  two-foot  thickness  of  wet-mix  shot- 
crete  for  a  distance  of  30  ft,  until  it  transitions  to  the 
natural  RCC  surface.  A  fixed  orifice  in  the  intake 
tower  regulates  discharge  to  a  maximum  of  60  cfs. 
This  discharge  rate  was  sized  to  drain  the  reservoir  in 
20  days  and  produce  minimal  flows  in  the  down¬ 
stream  channel.  An  intake  tower,  attached  to  the 
upstream  face  of  the  dam,  provides  inlet  control  for 
increasing  heights  of  sediment  deposition.  The  tower, 
a  typical  U.S.  Soil  Conservation  Service  design,  con¬ 
sists  of  a  double  weir  overflow  at  the  top  and  portal 
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intakes  at  5-ft  intervals  along  two  sides  of  the  tower. 
These  portals,  sealed  as  sediment  from  periodic 
impoundment,  accumulates  against  the  structure.  The 
structure  is  designed  to  be  an  unmanned  project. 
Discharged  water  flows  through  a  48-in.  outlet  pipe 
cast  in  the  RCC  dam  and  training  wall  into  an  impact 
basin.  Subsequent  low  velocity  flows  are  channelized 
and  eventually  discharged  into  the  natural  channel. 

d.  Floodway  channel. 

(1)  A  natural  water  course  below  the  dam, 
incised  into  the  canyon,  channels  water  flow  until  it 
reaches  the  city  limits  where  the  natural  channel 
widens  out  at  west  7th  Avenue  and  Vancouver.  At 
that  point  the  channel  improvements  consist  of  a 
concrete  intake  structure  with  trash  racks  and  an 
earthen  dike  to  funnel  flows  of  400  cfs  into  a  buried 
conduit  consisting  of  a  78-in.  reinforced  concrete 
pipe.  The  conduit  proceeds  east  on  West  7th  Avenue 
then  north  on  Rainier  Street  to  the  golf  course  where 
it  flows  out  from  a  concrete  stilling  well  and  follows 
a  natural  drainage  path  through  the  golf  course. 

From  Canal  Drive,  which  borders  the  golf  course  on 
the  north  side,  the  water  flows  through  a  6-  x  8-ft 
concrete  box  culvert  with  a  capacity  of  620  cfs  under 
Canal  Drive  to  the  Burlington  Northern  Railroad  fill 
where  a  78-in.  diameter  steel  culvert  was  jacked 
through  the  fill  to  be  able  to  pass  flows  up  to  620  cfs 
with  3-ft  of  freeboard  on  the  railroad  fill. 

(2)  Downstream  of  the  railroad  fill,  a  200-ft 
floodway  dike  was  constructed  to  elevation  383.5 
between  Highway  395  fill  to  the  high  ground  near  the 
Union  Pacific  Railroad.  An  opening  was  left  in  the 
dike  to  allow  train  traffic  to  continue,  with  a  stockpile 
of  material  near  by  to  fill  in  the  opening  when  flows 
exceed.  The  lower  Zintel  Canyon  water  course,  also 
know  as  Tweedt  Canyon  Drain,  is  a  combination  of 
natural  flow  channels,  low  bridge  crossings,  and 
culverts  crossings  under  embankments.  Depending 
upon  flow  magnitude,  water  will  either  flow  com¬ 
pletely  through  the  area  in  a  series  of  existing 
channels  and  culverts  or  escape  the  watercourse  and 
proceed  to  the  east  of  Highway  395  overpass  to  the 
Columbia  River. 

e.  Construction  operations  (Photos  C-1  through 
C-7). 

(1)  Crushed  basalt  rock  (140,000  tons)  for  the 
RCC  was  produced  from  a  quarry  located  only  a  few 
hundred  feet  upstream  of  the  dam  right  abutment. 


The  quarry  area  was  developed  using  a  dozer  (Cat 
D9C)  and  ripper.  A  crushing  operation  was  set  up, 
and  consisted  of  a  primary  jaw  crusher,  an  impact 
crusher,  and  two  roller  crushers.  The  RCC  mix 
required  29  to  32  percent  of  each  rock  product, 

2.5-in.  rock,  3/4-in.  rock,  fine  aggregate,  and  approxi¬ 
mately  6  percent  added  silt.  Approximately  50  per¬ 
cent  of  the  total  required  aggregate  was  produced 
prior  to  the  start  of  RCC  placement. 

(2)  Design  parameters  require  the  RCC  to  attain 
a  minimum  compressive  strength  of  1,400  psi  at  one 
year  of  age.  Static  stability  requires  cohesion  values 
of  35  psi  at  the  base  of  the  structure,  and  lesser  val¬ 
ues  in  the  upper  regions  of  the  dam.  Subsequent 
dynamic  analyses  determined  that  lift  joints  also  had 
to  attain  cohesion  values  of  50  psi  in  the  upper 
regions  of  the  structure.  It  was  determined  that  the 
specified  construction  system  had  to  provide  joint 
quality  that  resulted  in  shear  cohesion  values  exceed¬ 
ing  50  psi.  The  resulting  mix  attains  a  1-year  com¬ 
pressive  strength  of  2,200  psi,  and  displays  laboratory 
cohesion  values  of  95  psi  and  150  psi  for  unbedded 
and  mortar-bedded  lift  joint  configurations,  respec¬ 
tively,  at  exposures  of  24  hr  at  70°F.  The  paste-to- 
mortar  ratio  is  approximately  0.50,  the  mortar  volume 
is  23  percent,  and  the  workability  level  is  approxi¬ 
mately  15  sec,  measured  using  the  modified  vebe 
apparatus. 

(3)  Great  economy  is  achieved  when  RCC  pro¬ 
duction  and  placement  can  proceed  uninterrupted  at  a 
consistent  production  rate.  Repeatedly  changing  mix 
designs  (e.g.,  for  upstream  and  downstream  richer 
RCC  zones)  creates  placing  problems,  and  limits 
equipment  selection.  Consequently,  only  a  single 
RCC  mixture  was  produced  for  Zintel  Canyon  Dam, 
so  that  continual  plant  changes  were  not  required. 

This  is  especially  beneficial  for  continuous  mixing 
operations,  since  there  is  usually  no  convenient 
method  of  instant  and  frequent  mix  changes.  Several 
other  mixes  were  used  on  the  project.  A  higher 
cement  content  mix,  with  an  air  entraining  admixture, 
was  used  for  the  top  two  ft  of  the  stilling  basin  slab, 
as  well  as  for  the  top  four  lifts  of  the  dam.  A  low 
cement  content  mix,  with  an  air  entraining  admixture, 
was  used  for  the  top  four  lifts  of  the  training  walls. 

(4)  Precast  panels  for  vertical  face  construction 
were  fabricated  in  a  commercial  precast  facility 
100  miles  from  the  project  and  then  trucked  to  the 
site.  The  panels,  4  ft  by  16  ft  in  width  and  4  in. 
thick,  were  keyed  along  the  horizontal  joints.  The 
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panels  were  anchored  into  the  RCC  with  8-ft  coil 
rods  (3/4-in.  diameter)  and  end  plates.  Panels  were 
used  for  the  vertical  faces  of  the  stilling  basin  training 
walls  and  for  the  above-grade  vertical  surfaces  of  the 
upstream  face  of  the  dam. 

(5)  Panels  were  placed  in  a  checkerboard  config¬ 
uration  so  that  intermediate  panels  were  supported  by 
previously  placed  and  anchored  panels.  This  elimi¬ 
nated  the  need  for  external  bracing.  The  checker¬ 
board  method  of  panel  installation  is  a  very 
economical  panel  system,  however,  tight  alignment 
tolerances  are  difficult  to  achieve.  The  specified 
alignment  tolerances  were  purposely  broad  so  that 
such  a  panel  installation  system  could  be  utilized. 

(6)  The  sloping  surfaces  were  to  be  a  free- 
formed  RCC  slope.  In  order  to  dress  these  slopes, 
the  free  slopes  had  to  be  trimmed  with  a  backhoe 
bucket  periodically.  This  produced  the  relatively 
uniform  appearance  of  the  slope,  and  removed  the 
uncompacted  RCC  on  the  slope. 

(7)  RCC  was  conveyed  from  the  plant  to  the 
placement  and  discharged  directly  into  front-end 
loader  (Cat  980)  buckets.  The  material  was  driven  to 
the  specific  placement  location  and  dumped  onto  the 
uncompacted  RCC  surface.  The  dozer  (Cat  D4) 
spread  the  material  in  14-in.  thick  layers.  Compac¬ 
tion  was  done  with  a  10- ton  double  drum  vibrator 
roller  (Ingersol  Rand  DA-50),  and  supplemented  with 
a  smaller  roller  (Ingersol  Rand  DH-22).  Edge  com¬ 
paction  was  done  with  a  rammer  (Wacker).  Since 
Zintel  Canyon  Dam  required  only  moderate  shear 
performance  at  the  lift  joints,  bedding  mortar  was 
applied  to  the  lift  joints  to  assure  shear  and  tensile 
strengths,  and  vehicle  transportation  on  the  surface 
was  allowed  to  reduce  project  costs.  This  arrange¬ 
ment  did  not  jeopardize  the  lift  joint  quality  and  still 
provided  significant  equipment  cost  savings. 

(8)  RCC  was  placed  in  lifts  12  in.  thick  and 
mortar  was  applied  to  each  lift  surface.  To  minimize 
the  impacts  of  mortar  application,  the  contractor 
formulated  a  system  to  pump  mortar  to  the  placement 
and  shoot  the  mortar  on  the  surface.  The  mortar  mix 
was  modified  with  “a  high  range  retarder”  to  produce 


phenomenal  extended  set  times  and  reasonable 
strength  performance.  This  process  proved  to  be  very 
effective  in  reducing  manpower  dedicated  to  mortar 
placement  and  provided  uniform  coverage  of  mortar. 
The  retarder  is  a  product  originally  developed  to 
delay  the  setting  of  concrete,  in  transit  mixers,  for 
extended  periods  of  time. 

(9)  Placement  began  in  the  key  trench,  with  a 
placement  of  16  lifts,  totaling  1,800  cu  yd.  The  RCC 
was  conveyed  to  the  placement  and  dropped  to  the 
rock  or  RCC  surface  by  elephant  trunk  followed  by 
dozer  spreading  and  compaction.  The  placement  area 
then  expanded  to  the  stilling  basin  slab,  with  12  lifts 
averaging  1,400  cu  yd.  RCC  was  conveyed  to  load¬ 
ers,  and  subsequently  transported  to  the  placement 
location.  Loaders  traveled  as  much  as  possible  on 
fresh  RCC  surfaces  rather  than  the  older  surfaces  that 
were  being  prepared  for  the  next  lift.  Upon  comple¬ 
tion  of  the  stilling  basin  slab,  the  placement  area 
narrowed  to  84  ft  and  continued  to  narrow  as  the 
dam’s  height  increased.  The  RCC  lifts  for  the  stilling 
basin  training  walls  were  placed  concurrently  with 
each  lift  of  the  dam  placements. 

(10)  Production  rates  averaged  50  cu  yd/hr  dur¬ 
ing  the  early  key  placements  and  the  upper  lifts  (in 
the  upper  section  of  the  dam).  Typical  production 
rates  of  200  to  225  cu  yd/hr  were  maintained  during 
placement  of  the  stilling  basin  and  main  dam  lifts. 

The  typical  placing  sequence  was:  1)  vacuum  accu¬ 
mulated  debris,  ponded  water,  and  segregated  aggre¬ 
gate;  2)  air  clean  the  surface;  3)  wet  the  surface; 

4)  apply  bedding  mortar;  and  5)  place  the  RCC. 

(11)  A  drilling  program  commenced  approxi¬ 
mately  6  months  after  completion  of  the  RCC  place¬ 
ments.  The  purpose  of  the  drilling  was  to  remove 
6-in.  diameter  cores  from  the  structure  and  the  foun¬ 
dation  to  evaluate  the  actual  engineering  properties  of 
the  RCC  and  the  foundation  rock.  This  testing  pro¬ 
vided  excellent  information  for  future  design  efforts 
using  RCC.  The  testing  showed  that  shear  cohesion 
of  the  RCC  lift  joints  more  than  doubled  with  the  use 
of  bedding  mortar  on  the  lift  surfaces  from  85  psi  for 
unbedded  lifts  to  205  psi  for  bedded  lift  joints.  The 
parent  RCC  containing  no  lift  joint,  tested  at  290  psi. 
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Figure  C-2.  Site  plan 
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Photo  C-6.  Constricted  placement  conditions  for  RCC  placement 
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